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We model both the photo-stationary state and dynamics of an illuminated, photo-sensitive, glassy 
liquid crystalline sheet. To illustrate the interplay between local tilt 6 of the sheet, effective incident 
intensity, curvature and dynamics, we adopt the simplest variation of local incident light intensity 
with angle, that is cos#. The tilt in the stationary state never overshoots the vertical, but maximum 
curvature could be seen in the middle of the sheet for intense light. In dynamics, overshoot and 
self-eclipsing arise, revealing how important moving fronts of light penetration are. Eclipsing is 
qualitatively as in the experiments of Ikeda and Yu (2003). 


I. INTRODUCTION 

Nematic networks change their shape when their ori¬ 
entational order is induced to change thermally or, if dye 
molecules are present, optically. Mechanical strains can 
be between several hundred % for elastomers and e.g. 
4-5% for glassy networks. We concentrate on the latter 
since they are strong (moduli ~ 10 9 — 10 10 Pa) and their 
directors immobile, eliminating further causes of mechan¬ 
ical response. The experiments of Ikeda and Yu [1] show 
two uniquely interesting new phenomena: 

(i) Large, optically-driven response in the direction of 
the polarization of light. Samples are polydomain so 
the only definers of directions are fco, the incident light’s 
wave vector, and E 0 , the incident light’s electric (and 
polarization) vector. At once, both the ease of deliv¬ 
ery of stimulus and the control of mechanics by optics 
were demonstrated. These aspects have been explored 
by many authors [2-7]. 

(ii) Curling, in the direction of Eq, occurred with large 
amplitude. Remarkably, curling continued so that much 
of the photo-glass sheet eclipsed itself; see Fig. 1. In 
shadow, one might expect curling to cease, or indeed re¬ 
verse since the initially lower side of the sheet is now 
uppermost and being irradiated. No explanation of this 
seemingly paradoxical phenomenon was advanced by [1] 
or by subsequent authors. Our theory here shows that 



FIG. 1. A nematic sheet of Ikeda and Yu bending in response 
to illumination from above and self-eclipsing as the deforma¬ 
tion develops. The right hand half is stuck to the support. 


the phenomena uncovered by these seminal experiments 
reveal much about the non-linear and dynamical pro¬ 
cesses behind nematic photo-solid absorption of light and 
mechanical response. We suggest further experiments. 


II. ABSORPTION, PHOTOMECHANICS AND 
BEND ACTUATION 

Dye molecules are linear in their ground ( trans , t) 
states and bent when excited ( cis , c) by photon absorp¬ 
tion. The number fraction of cis, n c = 1 — n t , increases 
by illumination, I(x), and decreases by thermal recovery 
to nt at a rate 1/r, where r is the c-lifetime, and where 
I(x) is the intensity (Poynting flux) at depth x into the 
sheet/cantilever. Thus 

Oti c I 

r—A — n t -n c = h c = -(al/l 0 + l)n c + al/l 0 , (1) 

Ot i m 

where I m is a material parameter l// m = Tr and a = 
Io/Im = T/ 0 /(l/r) [8]. The constant T subsumes an 
absorption cross section per chromophore and a quantum 
efficiency, while the reduced time t/r derivative d/d(t/r ) 
is denoted by • . We neglect c-absorption, background 
absorption and scattering in order to simply establish 
the qualitative aspects of the Ikeda and Yu phenomenon, 
which we succeed in doing. Thus such additional sources 
of absorption are seemingly not central, and our model 
evidently has the essence of this mysterious effect. 

The parameter a measures the ratio of the forward rate 
T/ 0 (using the surface light intensity) to the back rate 
1/r. Large a implies strong perturbation from n c = 0. 
Small a is the Beer limit where n c ~ 0 and absorption is 
by a dye population little perturbed from the dark state. 
We show that the Ikeda and Yu experiments reveal non¬ 
linearity (a > 1) is vital. The photo-stationary state, 
h c = 0, gives 

al 

1+aT 


n, 


( 2 ) 



2 


where X(x) = I(x)/Iq is an intensity at depth x reduced 
by the intensity of light just having entered (at x = 0). 

Intensity is reduced with depth by the photon absorp¬ 
tion in eq. (1) above, leading to nt —> n c 


dl_ _ _rh , 
dx~ d 


( 3 ) 


compared with the decay rate 1/r - the penetration is 
linear and very deep, X ~ 1 — x/(ad) rather than expo¬ 
nential, X = e~ x / d , which accounts for substantial me¬ 
chanical response [8, 10, 11] when one would otherwise 
expect little response due to only a thin skin d -C w be¬ 
ing Beer-penetrated. Dynamics will turn out to be strong 
evidence for non-linear effects. 


where the Beer Length d subsumes cross sections, num¬ 
ber densities etc., and absorption depends on the number 
fraction n t of absorbers. When n t = 1, (q < 1), the Beer 
limit I[x) = Ioe~ x / d is obtained. We require finite con¬ 
version to get a mechanical response at all and to obtain 
the observed dynamics. Hence (3) must be solved in the 
non-linear limit of n c (x) ^ 0, that is, n t is a function of I 
itself [8], either statically, eq. (2), or dynamically, eq. (1). 

Creation of cis isomers lowers order and gives a photo¬ 
contraction along E, that is a strain e p = —Cn c in 
its simplest form, with C a dimensionless scaling. For 
e p ~ —0.04 and n c ~ 0.8 (say), then C ~ 1/20. An e p (x) 
varying with depth gives curvature 1 /R as the solid aims 
to reduce the elastic cost of deviating from its new, nat¬ 
ural local shape. The effective strain is 

e(x) = ^+K-e p . (4) 

The longitudinal stress at a depth x in the sheet is 
a modulus times this strain. Integrating the stress and 
the moment of the stress through the thickness, w, of 
the sheet to get the force and the torque, setting both 
these to zero, and cancelling the modulus yields the two 
equations [8] 

pW pW 

0 = / [— +K+ Cn c (x)]dx = / — +K + Cn c (x)]xdx. 

Jo R J o R 
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Eliminating between these two equations for w/R one 
obtains 


w 12 C f w fw \ . . , 

R = l^h l2“V n ' Wdl - (6) 

These equations must hold generally, even in a dynam¬ 
ically evolving system, in a limit where inertia can be 
ignored, for instance in the creeping motion seen by [1]. 


Taking the solution (7) for X into n c = aX/(l + aX) 
and putting n c (I(x)) into eq. (6), then integrating using 
a variable change d{x/d) = —(1 /X + a)dX, and finally 
eliminating for 1/R gives [4, 8]: 

^-(l-X.)(l-^)-f(l-^)] 

2 (8) 

where the dimensionless combination D = 12C(—)“ sets 
a scale to the reduced curvature w/R. (For C = 1/20 
and w/d = 3 of our illustration, D = 1/15.) Recall the 
definition of the reduced intensity, X , after eq. (2). The 
solutions of (7) for x = w are injected in this curvature 
expression. The curvature 1/R increases as incident in¬ 
tensity, measured by a, increases, but eventually must 
decrease again for very high a - penetration is deep and 
X(w) ~ 1 — w/(ad) —» 1. The cis fraction n c saturates to 
a high value (dye is depleted) and the consequent small 
variation of photo-strain e p (x) with depth cannot induce 
bending. 




III. PHOTO-STATIONARY DYE 
POPULATIONS AND MECHANICAL RESPONSE 

Using the stationary population (2) in the form rit = 
1 — n c = 1/(1 + al) in (3) for X gives dX/d(x/d) = 
—X/(l + aX). Integration gives [8, 9] 

ln(X(x)) + a{X{x) — 1) = —x/d. (7) 

The solution for X(x) is in terms of the Lambert-W 
function (or ProductLog function), W{c), which satis¬ 
fies c = W(c)e w (°\ Thus X(x) = V(ae a ~ x ^ d ). For 
large a - intense light giving a high forward (4 c rate 


FIG. 2. Curvature as a function of incident light intensity a 
for w/d = 3, with the optimal intensity, a m at <C>, for this 
thickness indicated. Two intensities a = 1,5 (x,o) are used 
illustratively below. 

Fig. 2 shows the non monotonic function of curvature 
against incident light intensity a: When incident light is 
very weak, the cis fraction is small thereby producing a 
small curvature. When incident light is very strong and 
therefore deeply penetrating, the cis fraction is close to 1 
and nearly uniform through the thickness, and thus again 
there is hardly any curvature. Curvature is maximized at 
an optimal intermediate reduced light intensity, a m say. 
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A. Photo-stationary shapes of sheets 

We are concerned with the interplay between bend as 
a function of incident intensity, but also bend making the 
sheet oblique to the incident light and thereby itself influ¬ 
encing the effective intensity and penetration of the light 
responsible for the bend. This connection and feedback 
determines the photo-stationary state and also the com¬ 
plex dynamics that we later examine. To establish and 
illustrate such qualitative effects, we adopt the simplest 
possible variation of intensity penetrating the upper sur¬ 
face of the sheet, namely / 0 cos 0 which simply expresses 
the dilution of the incident flux Jo to an effective flux 
by obliquity. This assumption will be deficient in detail 
since (a) there are complicated Fresnel coefficients gov¬ 
erning the wave amplitude refracted into the medium, 
and (b) the wave entering is not completely refracted to 
be along the normal to the sheet (though it is nearly so 
in the case of strong absorption where Snell’s law takes 
an extreme form). In consequence, the electric vector 
will not be along the local beam direction, see Fig. 3(b), 
and its effect on inducing bend will have another angu¬ 
lar factor. Consequently, though our analysis reproduces 
the qualitative features of self-eclipsing and dynamics de¬ 
pendent on trans-depletion front penetration, the precise 
shapes of sheets to be compared with future experiments 
may depend on these refinements of angular dependence. 



FIG. 3. A cantilever illuminated from above by light entering 
the solid with intensity Iq. (a) Initially flat, (b) curled up so 
that the tangent t(s) at s makes an angle 9(s) with y , and 
light is incident at angle 9(s) to the cantilever’s normal. 


At an arc distance s along the sheet, the tangent and 
normal are rotated through an angle 0(s), see Fig. 3(b). 
It is the local x direction (the thickness direction) that 
enters the attenuation equation (7) and the effective in¬ 
tensity is aocos0, where op = Io/I m is the effective in¬ 
tensity were the beam to strike normally. Hence 0(s) 
enters the solution I to be injected into eq. (8) for the lo¬ 
cal curvature, itself having explicit 6 dependence through 
the ao cos 0 factor appearing. 

Eq. (8) can be written in a superficially simpler form 
as 

d0/d s = = — acos0(s), (9) 

R w 

where there is also 0-dependence in a: 


d — Op [ ^ '-^w ( 1 -E- w ) ( 1 



Q!o COS 9 
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both explicitly, and buried in X w . This differential equa¬ 
tion can be integrated to give 0(s), and then a second in¬ 
tegration gives the photo-stationary shape (x(s),y(s)) of 
the sheet bending in the x — y plane. See Fig. 4 for shapes 
corresponding to two effective intensities ao (that would 
be falling on flat sheets) that are less than or greater than 
a m ; the two cases are qualitatively different. 



FIG. 4. Photo-stationary shapes for ao — 1 and ao = 5. 
The positions of maximum curvature correspond to angles, 
and therefore arc positions, marked on Fig. 2 (x,0). The 
reduction of curvature by w/D is used here for lengths too. 
Sheets of reduced length L/(w/D ) = 14 are shown. For a 
given ao and w/d, these are master curves; sheets of smaller 
L will terminate at intermediate places along the curve. 


Consider incident light with ao = 1 < Qf m , that is 
smaller than the a for optimal curvature for this w/d. 
The effective intensity incident locally on the sheet, 
a = aocos0, decreases with 9. Therefore for a 0 = 1, 
the maximum curvature obtains at s = 0, as shown in 
Fig. 4, where a(s) is greatest, namely ao- However inci¬ 
dent light with ao = 5 is more intense than optimal: As 
9 increases, a decreases down to a m (if the total length L 
is large enough) and curvature increases to a maximum. 
If 9 continues to increase, then a(s) and thus curvature 
decreases to 0. Thus the maximum curvature for ao = 5 
obtains at an intermediate s in the sheet; see Fig. 4. 

To give a simple insight into the effect of curvature 
being arc-position dependent, one can look for solu¬ 
tions of eq. (9) for curvature where we take the ex¬ 
plicit, leading 0-dependence via cos0, and ignore the 
further 0-dependence in a which we now set to be con¬ 
stant. Reducing arc lengths by w/(aD), that is s = 
uw/(aD), then d0/du = cos 9(u) integrates simply to 
sin0(u) = tanh(u). Recognising that (x(u),y(u)) = 
w/(aD ) J o “d« , (cos0(i( , ),sin0(?/)) and using cos 9{u) = 


2 
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dO/du yields the parametric forms. 

/ \ W /!/ x W . -1 / W XX 

xm) = •— —u(u) = — sin (tanh(u)) 
au aV 

w 

y(u) = — ln(cosh(it)) (11) 

and eliminating u gives 
w 

y(x) = --In (cos \x/(w/aD))) (12) 

aD 

which is qualitatively of the form in Fig. 4. 


B. The marginal effect of gravity 

To rule out any role of gravity in determining static 
shapes, and later the dynamics of a glass sheet or can¬ 
tilever, we calculate an extreme bound from Fig. 3(a). 
Considering such a sheet of width W, we can calculate 
the curvature close to the clamped end induced by grav¬ 
ity, the torque being estimated as if the sheet were flat 
(i.e. estimating the most critical location for curvature 
and over-estimating the torque). Equating elastic and 
gravitational torques, 

EWw 3 pwWgL 2 w pg 6L 2 _ 6 L 2 
12 R 2 ^R~lE^ur = lw !l '‘ 


over-bend, and this is further complicated by (a) the flux 
driving the depletion front varies with 9(s), which results 
from the accumulation of curvature 89/8s from s' = 0 
to s, and (b) development of angles 9(s,t) > 7r/2 giving 
eclipsing, the incident light being blocked from falling on 
the sheet at arc positions before that s where 9 = 7r/2. 
Such sections are in the dark, their n c {x) fraction of cis 
recovers, and they lose their curvature. The sections of 
sheet that are doing the eclipsing necessarily have their 
under sides now exposed to the light and their curva¬ 
ture is reduced and eventually reversed. A reversal can 
lead to double eclipsing. We now explore this complex 
dynamics. 

The cis fraction at any time can be taken from (3), 
solving for n t and replacing it by 1 — n c . Thus n c = 
1 + (d/I) 8I/8x = 1 — ddA/dx where the absorption 
A = — In (I). When the above n c is injected into eq. (6), 
the 1 term in n c gives a vanishing integral. The gradient 
term either integrates trivially against w/2, or integrates 
by parts against —x to give overall 


w 89 _ w/D 
Dds = R 


’-A(w,t) — [ A(x,t) 
Jo 


da? 

T‘ 


(14) 


We have used A(0, t) = 0 since 1(0, t) = 1. 

The coupled, non-linear pair of PDEs (1) and (3) for 
n c (= —ht) and I’ can be reduced to a single temporal 
quadrature for A [14] at each s: 


where l = E/pg is a characteristic length emerging from 
matching elasticity (E - the Young’s modulus) with grav¬ 
ity (SO P ~ acceleration of gravity and density of the photo¬ 
glass). For photo-glasses l ~ 10 5 m. 

If gravity were to change the beam’s curvature by 
that of a quarter circle of radius R , thus L = ttR/2, 
then 1/R = 7r/2 L in the above yields a length L g = 
(irlw 2 /12) 1 / 3 before gravitational effects compete with 
elastic. Ikeda and Yu [1] have w = 7pm and L ~ 3mm 
(half their total sample length since they clamped in the 
middle and not one end), whence L g ~ 10~ 2 m is comfort¬ 
ably larger than their L. Mol et al [12] have w ~ 10/im 
- 40/rrn, and their L = lCF 2 m is comfortably within the 
(extremely low) estimates of L g = 2 — 4 x 10~ 2 m. We 
henceforth ignore gravitational effects. 

IV. DYNAMICAL PHOTOMECHANICAL 
RESPONSE AND ECLIPSING 

The varying cis-fraction of dye with depth and time, 
n c (x, t), arises from penetration of a trans depletion front. 
Sometimes this front is loosely called a bleaching wave 
or front, since the converted dye is not an effective ab¬ 
sorber of the original colour of light when in the cis state. 
However this bleaching is reversible (not chemical) and 
is otherwise referred to as saturated absorption [13]. The 
resulting curvature (6) determined by n c is thus a func¬ 
tion of time, non monotonic since 1/R depends on the 
spatial variation of n c (x). Thus we have time dependent 


A = x/d — A + a 0 cos9(s,t)(e A — 1), (15) 

with A(0, f) = 0 and A(x, 0) = x/d. The latter is Beer’s 
law of exponential decay I = e~ x / d of light in an as-yet 
undepleted dye population. This 2nd condition needs 
careful reexamination after eclipsing when sheets start 
to be irradiated from the back face. 

We solve (15) (see [14]) for depletion front solutions 
and inject them into (14) for (w/D)/R(s,t). The form of 
the dynamics, especially eclipsing, depends critically on 
the sheet length L and the reduced intensity a. Longer 
L means that higher angles 9(s,t) = ds'89/ds' can 

be accumulated and 9 —>• 7r/2 is more achievable. When 
the tip of the sheet approaches 7t/2, it can be convected 
over to 9 > 7t/ 2 by continuing light penetration at s < L 
since curvature increases with a further increasing gra¬ 
dient of n c (x) deeper through the sheet thickness. This 
non-linearity in response leads to eclipsing before later 
recovery. 

Fig. 5(a) shows the dynamics for ao = 5. The corre¬ 
sponding photo-stationary state in Fig. 4 is the dashed 
line in Fig. 5(a). The sheet starts bending from an ini¬ 
tially flat shape. One already sees overshoot at t/r = 0.1. 
The overshot section is now illuminated on what was the 
back face, while at least part of the sheet with 9(s) < n/2 
is eclipsed. The maximum curvature is more or less 
in the middle of the sheet where it was most strongly 
bent before being eclipsed. After overshoot, parts with 
9(s) > 7t/ 2 bend backward because of the reversed il¬ 
lumination, as can be seen in t/r = 0.3, while eclipsed 
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FIG. 5. Bend curve sequences in t/r of 0.01, 0.1, 0.3, 4, oo 
for (a) ao = 5 and (b) q 0 = 1 (where the t/r = 4 curve sits 
on top of that for t/r = oo). Reduced sheet length is again 
L/(w/D) = 14, and w/d = 3. 


parts unbend exponentially in time since they are in the 
dark. Approach to the stationary state takes a long time, 
in terms of the fundamental time scale r, because of the 
complex sequence of overshoots. 


Fig. 5(b) has = 1, which is slightly less than the 
maximal value a m for w/d = 3. Now strongest bending 
is closer to the fixed end, as expected from Fig. 4 for 
this ao. Overshoot is not so extreme and the approach 
to stationarity is much quicker. For smaller «o, for in¬ 
stance 0.5 for this length L, overshoot still occurs, but 
only slightly, whereas for ao = 0.1 with this L, it is lost. 


V. CONCLUSIONS 

We have demonstrated that because the effective il¬ 
lumination is controlled by orientation, and at the same 
time drives orientation since it induces bend, the bending 
of photo-responsive sheets is subtle. To explore the qual¬ 
itative behaviour that should arise, both in statics and in 
dynamics, a simple geometrically-inspired dependence of 
light penetration on angle of incidence is adopted. The 
other essential physical driver of this photo-mechanics is 
that conversion of the isomerising guest dye molecules to 
their excited state has to be considerable (i) to perturb 
the local structure of the glass and induce mechanical 
change, and (ii) to allow deep penetration (via a front of 
depletion of the ground state species) so that bend can 
actually occur. We then find: 

• Photo-stationary shapes arising are qualitatively 
different according to whether the intensity of il¬ 
lumination normally incident is more or less than 
a characteristic value that depends on both a ma¬ 
terial constant and on the thickness of the sheet 
reduced by the Beer length for adsorption. These 
stationary states cannot be self-eclipsing. 

• Seminal experiments on nematic glass bend re¬ 
sponse did display self-eclipsing as response pro¬ 
ceeded. By analysing the dynamical evolution of 
the mechanics, we show that this must have arisen 
as a transient effect and, accordingly, must be im¬ 
possible to explain with just linear, Lambert-Beer, 
light absorption. We further predict that the route 
to the final photo-stationary state must show par¬ 
tial unbending via back bend, and could display 
multiple eclipses. 

The two regimes both suggest further experiments 
which should also give insight into how complex the true 
angular dependence of absorption is. For example the 
experiments of Ikeda and Yu [1] could be illuminated 
for longer, such that appreciable backbend (in addition 
to eclipsing) is observed. Eclipsing itself depends quite 
crucially on the length of the photo-responsive sheets, 
therefore repeating the experiments of Ikeda and Yu with 
different sample lengths should also reveal a complicated 
behaviour. Perhaps most simple of all would be to coat 
the back face of the sheet with a reflective coating, thus 
removing effects from illuminating the back surface and 
allowing one to focus entirely on front-surface illumina¬ 
tion and eclipsing. 
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